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A weak measurement consists in coupling a system to a probe in such a way that constructive 
interference generates a large output. So far, only the average output of the probe and its variance 
were studied. Here, the characteristic function for the moments of the output is provided. The 
outputs considered are not limited to the eigenstates of the pointer or of its conjugate variable, so 
that the results apply to any observable o of the probe. Furthermore, a family of well behaved 
complex quantities, the normal weak values, is introduced, in terms of which the statistics of the 
weak measurement can be described. It is shown that, within a good approximation, the whole 
statistics of weak measurement is described by a complex parameter, the weak value, and a real 
one. 



I. INTRODUCTION 



A weak measurement is but a peculiar interference ex- 
periment: a quantum system is initially prepared in a 
superposition of eigenstates of an observable A; the sys- 
tem interacts with a second quantum system, the probe, 
prepared in a suitable coherent state, through a weak 
coupling proportional to A; as a result, the probe evolves 
as if subject in parallel to different interactions; the dif- 
ferent branches of the probe's wavefunction interfere, as 
revealed by the observation of the pointer variable; as 
in an ordinary interference experiment the fringes disap- 
pear if averaged over the whole screen, the same hap- 
pens with a weak measurement; in order to fix a point 
of the 'screen' it is necessary to make a postselection 
of the system. Weak measurement (without postselec- 
tion) was introduced in the pioneering paper of Arthurs 
and Kelly pLi , where it allowed to overcome the quantum 
mechanical limitation to the joint measurement of non- 
commuting observables. The coherence of the detectors 
manifests in the variance of the outputs [2], resulting in 
an uncertainty relation ApAq > H, twice as large |3j as 
the Kennard limit [U [S]. In another pioneering paper 
[S], Aharonov, Albert, and Vaidman studied the aver- 
age output of a weak measurement of a single observable 
with post-selection [7], showing that it can be dramat- 
ically larger than the maximum eigenvalue. This effect 
was soon attributed to the coherence of the probe [S]. 
Interestingly, when unpolarized spins are injected into 
a series of three detectors, even though these are ini- 
tially prepared with no coherence in the readout basis, 
coherence can be induced in the middle detector by the 
interaction with each spin, while the first and third de- 
tector provide pre- and post-selection. This has observ- 
able consequences if the rate at which particles are fired 
into the detectors is higher than the decoherence rate of 
the latters [9^ . Several experimental realizations of weak 
measurement have been performed by now [TD^EO] . Weak 
measurement has been contemplated theoretically for ap- 
plications in solid state devices [21-27 . On a theoretical 
ground, Ref. [28] extended the results of [6] to any initial 
state of the probe; also, the average and the variance of 
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FIG. 1: A schematic view of the measurement with pre- and 
post-selection, the horizontal direction representing increasing 
time. A system is prepared in pi, interacts with a probe, then 
is projected into 15*) by a measurement. After the interaction, 
the observable o of the probe is measured. A random number 
generator produces a number 0<a;<l. Ifa;> w{S), an 
arbitrarily chosen probability depending on S, then the record 
o is discarded. 

both p, the pointer variable, and of g, its conjugated vari- 
able were considered [5S]; in Ref. [5D], the general case 
of a non-instantaneous interaction was studied, showing 
that the dynamical phase of the detector influences the 
outcome. Ref. [5] associated a complex quantity, the 
weak value A™, to the weak measurement. Assuming a 
Gaussian distribution for the probe, the real part of the 
weak value equals the average output. For spin 1/2, this 
value can easily exceed h/2 even for random pre- and 
postselected states [3T]. The real part of A™ describes 
well the average output of the weak measurement as far 
as the pre- and postselected states are not almost orthog- 
onal. When they are nearly orthogonal, the output can 
indeed reach a large value, but the perturbative expan- 
sion used in |6j fails, for a fixed coupling strenght [37], 
since A™ diverges, while the observed output has a finite 
maximum just before the orthogonal configuration, then 
drops abruptly [301133]. In simple words, weak measure- 
ment with postselection relies on a conditional probabil- 
ity which, by Bayes' rule [32], is V{o\f) = V{o, f)/V{f), 
the ratio of the joint probability V{o,f) of observing o 
for the probe and making a final postselection /, over the 
marginal probability V{f) of making a successful postse- 
lection. Both probabilities can be expanded in powers of 
the interaction strength. For nearly orthogonal pre- and 
postselected state, however, the leading term in the de- 
nominator vanishes faster than the corresponding term in 
the numerator for the average JdooVio, f)/'P{f). This 
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artificial divergence should be treated by expanding both 
numerator and denominator up to the first non vanishing 
terms. This apparently straightforward prescription was 
missed so far by the community working with weak mea- 
surement, with the exception of [33 , where the correct 
expansion, however, was made only for the denominator. 

Furthermore, while the experiments provide the full 
probability distribution V{o\f), the theory generally con- 
siders only average values, and at times variances [^151] . 
Here, we provide the full counting statistics, i.e., the 
characteristic function for all the moments of any ob- 
servable 6 of the probe, conditioned on the postselection. 
We give a correct expansion working for any pre- and 
postselected state, and a somewhat simpler interpolat- 
ing formula working satisfactorily over the whole range 
of postselection. 



II. MEASUREMENT FOLLOWED BY 
POSTSELECTION 

We start by describing an ideal, instantaneous (von 
Neumann) measurement with postselection. Let us con- 
sider a quantum system prepared, at time ti—e, in a state 
PsiU—s) — Pi (preselection), with e an infinitesimal time. 
The system interacts, at time t^, with another quantum 
system, the probe, through Hint — ^hXS{t~ti)qA, where 
A is an operator on the system's Hilbert space, q on the 
probe's. The probe is prepared in a state pp{ti — e) = po- 
At time tf = ti + e a projective measurement of the ob- 
servable o is made on the probe. The procedure is partic- 
ularly interesting, when d = p is conjugate to q, [q,p\ = i, 
since then it is the observable of the probe that carries 
information about the measured quantity A. At time tf, 

a projective measurement of an observable 5/ of the sys- 
tem is made, giving an output S and leaving the system 
in the state \Sf: S). Then, given S, one keeps the record 
o according to an arbitrarily chosen probability w{S). 
This leaves the system in the postselected mixed state 
Pf = j:s^iS)\Sf: S}{Sf: S\/W, with W = Es^i^)- 
The procedure detailed above describes a measurement 
with pre- and post-selection and it is sketched in Fig. [T] 
Since the interaction is instantaneous, the problem can 
be treated analytically. The joint state for the probe and 
the system after the interaction is 



^5 = e^^^V^ ® Poe-'^^^ 



(1) 



or, in the basis of eigenstates of p and A, 

p{p,a,p',a';tf) = po{p- Xa,p' ~ \a')p,{a,a'). (2) 

The conditional density matrix for the probe given an 
initial preparation on poi a system preselected in pi and 
successfully postselected in p/ at time i/, follows readily 
after applying Born's rule to Eq. ([T]) 

Pp\o,^J = Trs {pf <S> 1 e'^^^p, ® poe-'^^^, (3) 



with Tr the trace, and the subscripts S referring to the 
system (in the following, the subscript P will refer to the 
probe). In the readout basis, 

P(p,p'|0,z,/) = ^po(p- Aa,p' - Xa)G{a,a), (4) 



where 0, i, f indicate the conditions (from now on we shall 
write only / and take 0, i as implicit conditions), and with 
the definition 

G{a,a') := j^{a\pf\a){a\p,\a), (5) 



while the normalization is 



W 



' Poe 



= j dp'^poip- \a,p- \a'){a'\pf\a){a\pi\a') 

a, a' 

dqpo{q,q)Z'-{Xq,Xq), (6) 



where Z"'{p,i') is the weak characteristic function, de- 
fined by 



(7) 



N is proportional to 7^(/), which is the probability of 
making a successful postselection in irrespectively of 
the value of p, or of the fact that p, or any other observ- 
able of the probe, was measured at all. 

The joint probability of making a successful post- 
selection and observing an eigenvalue o of 6 for the probe 
is then 



V{o,f) ^WTrs.piPf ^^oe'^^'^p^® PQi 



-iXAq 



}, (8) 



withlio = |o)(o|. The conditional probability is 7^(0!/) — 

r{oj)/p{f). 

Eq. (|3| can be conveniently rewritten in the q- 
representation 



piq,q'\f) ^N~'po{q,q')Z^iXq,Xq'). 



(9) 



Equation ([9| holds also when q does not have a conju- 
gate variable. In this general case q may have a discrete 
spectrum, and then the integrals are to be interpreted as 
Lebesgue-Stieltjes ones. The exact conditional probabil- 
ity for any observable o is then 



7^(01/) = - j dqdq\q'\no\q)Z^{Xq,Xq')po{q,q'). (10) 

The characteristic function is obtained by the substitu- 
tion exp [ixo] — > Ho in Eq. (10), 



^(xl/) = ^ / dqdq'{q'\e'nq)Z''{M,W)po{q.q')- 



(11) 
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Equation Q shows the pecuUarity that the statistics 
of V{p\f) = p{p,p\f), is influenced by the off-diagonal 
elements of poip^p'), so that interference effects may be 
relevant. By contrast, Eq. ^ reveals that no interfer- 
ence effects arise if one chooses to observe q rather than 
p, so that this case is of lesser relevance. Interference 
terms are appreciable when XSa <^ Sp, where Sa is the 
spacing between eigenvalues of A and Sp is the coher- 
ence scale in p, i.e., the value over which the off-diagonal 
elements p{p + Sp,p — Sp) die out. A measurement per- 
formed in this regime is called weak measurement. The 
state of the system after a weak measurement is almost 
unaffected, but may acquire a small perpendicular com- 
ponent to its initial state. In the rare cases when the 
postselection is made on an almost perpendicular state, 
so that Tisipf Pi} — 0, the average value of p can be fi- 
nite, due to the constructive interference represented by 
the contribution of the off-diagonal elements. Since for 
a strong measurement there is a one-to-one correspon- 
dence between the final state of the detector and that 
of the system, so that A = p/X is associated to the sys- 
tem, if one makes stubbornly the same position, then the 
weak measurement of A can yield a large average output. 
This is a purely quantum effect, since if it was due simply 
to a large classical uncertainty in the initial distribution 
Po(j>,p), then each individual measurement may give a 
large value, but this would be washed out by averaging. 



is a projection operator (over the subspace spanned by 
the eigenvectors of A with non-zero eigenvalue). Spin 
1/2 and spin 1 (with n — 2), and yes- no observables, 
i.e., projection operators (with n = 1), are important in- 
stances. In this second case, a!;^„ vanishes as fast as a^g 
for NOPPS if either or is a mixture of states belong- 
ing to the subspace defined by A" (which may coincide 
with the whole vector space if A" = 1); consequently. 



<^o,o/ao,kn 1 and a 



j,kn 



for such states. However, 
i 7^ for NOPPS that do not fah under the previous 
case, so that the conditional probability is not indeter- 
minate. 

For brevity we put ajg = aj and in particular 
a™ = a™. The normal weak values appear naturally 
in the expansion of N and p{p,p'). For historical rea- 
sons, the canonical weak value used in the literature 
is A'^ = af/a^, which diverges for NOPPS. A recent 
paper 
A"! 



33] introduced a family of diverging quantities 
a^}./aQ in terms of which to describe the weak 
measurement. As this approach may obscure the rel- 
atively simple problem at hand, we shall use the well 
behaved normal weak values. 

The normal weak values (and the canonical weak 
value) arise naturally whenever the initial density ma- 
trix for the detector \p{q^q')\ is peaked around a value 
{q*;q*) such that \q* <C 1. If this is not the case, then 
one should define 



III. NORMAL WEAK VALUES 

The results presented so far are exact, but not always 
susceptible of an easy evaluation. In the following, we 
shall make an expansion in the interaction A, assumed to 
be small compared to some appropriate scales that will 
arise in due course. A set of quantities emerging naturally 
from the expansion are the normal weak values, defined 
as 



„ ^ {^^y+^z-'{^Ji,v) 



i/=0 



TrsiA^PfA'^p,}. (12a) 



We notice that a^^ — ctjl, so that ajj are real. For 
nearly orthogonal pre- and postselected state (NOPPS), 



a^ Q — >■ and uJq — > 0, in such a way that a^g/a^Q — ?> 



for j > 1; also, a^'^. -/^ for j,k > 1. There are two 
exceptions to this behaviour: 

1) Either pi or pf is a mixture of degenerate eigenstates 
of A. This case bears no interest, since one is basically 
repeating the same measurement, the second time in the 
weak regime, the first or the third in the strong one. 
Thus, it is impossible to observe a postselected state 
which is orthogonal to the preselected one, which refiects 
in all a™ J, being zero, so that the conditional probability 
takes the indeterminate form ^. 

2) Neither pi nor is a mixture of degenerate eigen- 
states of A, but A is n-idempotent, = A for some 
n > 1 (up to a factor |Mj)- This implies that = A^" 



„ ^ (,a).-+fcZ"'(^,z/) 



u—Xq* 



--Tr s {A^Pf A'' e'^'>'^p,e 



■i\q* A 



}■ 



(12b) 



These values are no longer specific to the system, but 
depend also on the coupling and the detector property 

q*- 

In order to express the characteristic function later on, 
we define a one-parameter family of weak values 



aife(2) := 



fl — Z 
V — — Z 



dfi^{—diyy 

r /i 7 aI^ izA izA^ 

=Tts{A-' pfA e p^e |, 
or, if Xq* is appreciable. 



(13a) 



dp^{—dvy 
=Trs{iV/i'=e*(^+^'*)^p,e*(^-^«*)^}. 



f_L—z-\-Xq'' 
ly— — z-\-Xq* 



(13b) 



IV. CHARACTERISTIC FUNCTION FOR q 

We proceed to evaluate the full counting statis- 
tics of V{q\f), i.e., its Fourier transform Zq{x) = 
J dqexp[ixq]'P{q\f). We need first to expand N, which 
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from Eq. Q (see also, e.g., [321) results 



with 



ri=0 



j=0 



(14) 



where the bar denotes average with po. Then, after ex- 
panding in Eq. ( 11 ), we have 



n=0 ■ j=0 



a A 



(15) 

In order to have an expansion working for NOPPS, we 
expand both the numerator and denominator to second 
order, which contains the non- vanishing quantities a^i 
(notice that Ref. [33J, when estimating the average, in- 
consistenly expands the numerator up to the first order 
and the denominator to the second order, so that the 
case of orthogonal pre- and postselected states has to be 
treated as a separate case, unnecessarily) 



Zq{x) =iV2"^|e'X9a^ -I- 2A9e*x«lm(a5") 

+ A2^V^(ar^i-Re(a^))l, (16) 



with exp [ix<i] the initial characteristic function and 
iV2 - ao + 2A9lm(a5") + A^^ [al^ - Re(a^)] . (17) 

After considering that higher order terms in A become 
important only for NOPPS, so that a™i ^ can 
use a simpler interpolation formula within a satisfactory 
approximation 



Zq[x) - ^2 e^^'ao + 2\qe-x<ilu,{a^) 



A2g2gzx9a5"i 



(18) 



+ 2AgIm«) + A2g2Q,^^. 



(19) 



However, if one is deriving the j-th moment (q-')/ from 
Zq, one should keep the term in the numerator of 



Eq. (16) whenever qi+^ — or gJ+i « Xqi^"^, since then 



the first order contribution vanishes or it is comparable 
to the second order one. 

In all the equations, one may divide both numerator 
and denominator by , so that the simpler interpolation 
formula is a function of a complex number, the canoni- 
cal weak value. A'" — Tr {p/Api}/Tr {p/pi}, and a real 
number B"^ = Tr {p/^lpiA}/ Tr {p/pi}. The inequality 
> 1^1" 1 2 jiolds, with the equality sign for pure pre- 
and postselected states. For NOPPS B"^ diverges faster 
than , dominating the expansion. When, as discussed 
above, q^+^ = or qi+^ Ag-'+2, an additional complex 
parameter should be used = Tr {p/A^pi}/ Tr {p/pi} 
in order to get a good interpolation for (q^) f. 
V. CHARACTERISTIC FUNCTION FOR THE 
READOUT 



A. Special case: Gaussian probe 



Let us say that the probe is described by an initial 
Wigner function 



WoiQ,P) 



1 



27rAPAg 



exp 



p2 



2AP2 



2AQ2 



(20) 



We recall that the coherence scale Sp in the readout basis 
is 6p = 1/2AQ. The exact expression for the probability 
of observing p is 



np\f) 



1 



2ttAP 



^exp 



iXq{a — a') 



A2AQ2 



f\2 



a') 



exp 



[p-X{a + a')/2f 
2AF2 



G{a,a'). 



(21) 



While one can expand the first exponential in the nu- 
merator in Xq, XAQ, provided these are small, it would 
be a grave error to expand the exponential containing p 
in A/AP, as it is appreciable only when p takes values 
order A. In order to calculate the moments of p, {p') f = 
J dpV{p\f) p^ one should change the integration variable 
top' ^p- X{a + a')/2, so that p> = [p' + X{a + a')/2]J, 
and expand the binomial. By following the above pre- 
scription, and keeping terms up to A^, we obtain, for 



even j , 

{p>)j ^pJ+^(^^^p 
and, for odd j, 



,J-2 



Re(a^)]/7V2, (22) 



{p')f =X [jy-iRe«) + XjW-nmia^ /N^, (23) 

The bar means average with respect to Pq (p) ■ After con- 
sidering that higher order terms in A become important 
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only for NOPPS, so that a'^i ^ a™, we can use within 
a satisfactory approximation 



= < + 2Xqlm{a\") + X'^q^a'^'^. 



(24) 



We may now divide both numerator and denominator by 
, in order to eUminate one parameter, as done in the 
previous section, obtaining for even j. 





{p>)f =pJ + 



)pj-2 [B"" +Re(C"')] 



2 1 + 2AqIm(A«') + X^q^ [B"' - ReCC"")] 

(25) 



and, for odd j, 



{p>)f ^Xj 



pJ-iRe(yl"') + Agpi-iIm(C"") 
1 + 2Xqlm{A^) + XW [B'^ - Re(C™)] ' 



(26) 



The characteristic fonction, defined as Zp{x\f) 



(27) 



Zpixlf) =eJXPZs{Xx) 
with g-x'APV2^ j^j^j 

Zsi^x) = ^e^^['^('^+"')/2+«>-'^')l-^("-"')'G(a,a') 

a, a' 

(28) 

In particular we have that the n-th cumulants (logarith- 
mic derivatives of Z) are of order A" for n > 3. After 
expanding in Xq, XAQ , 



e^xp < / Ay 



A2fl2 



Ax 
2 



2AqIm (^aj" 
- Re ( a^' 



Ax 
2 

Ax 



(29) 



and, up to second order terms in A, 



Zp{x) 



+ X' 



a^ + X [ixRe(a5") + 2qlm(a^)] 



4 

,.2 



Re(a^) 



(30) 



One should take Eqs. (22 1, (23), and ( |30| with a grain of 
salt: in Eq s. ^ and^23|, we have neglected terms of 
order (^)AV"~''(a + a')''/2'' compared to (2)AV"^(a + 
a')'^/4: for k > 2. Without loss of generality, if A has 
a limited spectrum, we can rescale A — > Aniax|a|, A — 
A/ max|a|, so that \a + a'\/2 < 1. For large enough n 
the approximation fails, since 



Cn„n-2 

n->-oo ("jp"~4 



Considering that (^)(p/AP)" '^x''' has a maximum for 
k 



'nx^, we have that the validity of Eqs. (22 1 and 
AP^/A^. The characteristic func- 



(231 extends to n* 



tion in Eq. ( 30 ) should be interpreted as giving the cor- 



rect moments up to n* . Since all the moments of p are 
necessary to reconstruct the probability distribution, one 
would not get the correct V{p\f) by Fourier transform- 
ing Z{x\f)- On the other hand, Eq. (29l provides a 



good approximation to the moments of all orders, even 
though the expansion in A is not done consistently. Once 



Eq. ( 29 ) is Fourier-transformed to the conditional prob- 
ability V{p\f), it captures the shifts p ^ p 



Xa induced 



by the interaction, contrary to Eq. (301 



B. General case 

The former particular case suggests the proper treat- 
ment of the expansion in A. We define the function 



Ro{P,p) ■.^Po{P + p/2,P-p/2) 



(31) 



Substituting in Eq. (|4| for p' — p, we have the conditional 
probability 

V{p\f ) = ^i?o(p - A^, A(a' - a))G{a, a'). (32) 



The correct way to proceed is to expand Rq in its sec- 
ond argument only. For the Gaussian distribution of the 
preceding subsection, Rq factors into the product of a 
function of {p + p')/2 and one of p — p' . This is not gen- 
erally the case. For the characteristic function, we have 

Zpixlf) = E / dpe'^^Ro (p - A^, A(a' - a)^ G{a, a') 



Ax 



iN~^\ e'XPa^ ['-Y) + '2Xe'^Pq Im ( i 2 



Ax 



+ X^e^xpqi 



Ax 



Re I 



Ax 
2 



(33) 



where the bar denotes quasi- averages, f{p)g{q) — 
J dpdq f{p)g{q)Wa{q,p)- In terms of properly defined av- 
erages, 

= ^ E (fc) V (P) (34) 



We notice that whenever /(p) = const, or g{q) — const. ^ 
the quasi-averages coincide with proper averages, so that 
there is no notation clash with the former section and 
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subsection. The moments, to order A^, are 



{pi) f =pi + X jy-iRe(a5") + [piq-p^q] Im(a5") 



(p^^-p^ q^) [ar,i - Re(a^)] 



(35) 



while the characteristic function can be approximated, 
with the same proviso as in the former subsection, as 



Zp{x\f) o,N^'\e^xPa^ 
+ X 



2 



q2e^xr>^^e^xp]al,+^xqe^xPlmia-) 



g2gixP + ^e»xp ) Re(a^) 



(36) 



As for the Gaussian case, the vahdity of Eq. (35) is 



VI. CHARACTERISTIC FUNCTION FOR ANY 
OBSERVABLE 



A. Small Xq* 



Finally, the results are easily extended to the case when 
one chooses to measure a generic observable 6 of the de- 
tector, instead of p or q. After expanding Eq.(ll), as- 



suming that Xq* ^ 1 (we recall that q* indicates the 
maximum of the distribution paiq, q)), we have 



with {X)o = Trp {Xpo}. Expanding up to second order 
terms. 



^i([q,e'>^°])oRe«) 



+ ({g,e'>^°})oIm«) 

^2{qe^^'q)oaZ, - ({g^, e^'^°})oRe(a^) 



(38) 



Putting 6 = g or 6 = p in Eq. (|38[), we recover Eqs. ( 15 ) 



and ( 36 1 . We notice that Eq. (p38| should be interpreted 



as generating formulas for the moments (o-' ) / , up to a 
critical n* that must be determined case by case accord- 
ing to 6 and po. 

As an example, we have that the expectation of 6 is 



-i[(7,6]Re(a5")-f {g,d}Im«) 



' 2 



2qdqal, - {e, 6}Re(«^) - 6]Im(a^) 

(39) 



If qdq 7^ one can neglect the terms in the second 
order, giving 



limited to j < n* . The estimation of n* depends 
on the actual form of the initial density matrix of the 
probe, and now it relies on comparing {'l)qp-'~^, 
and (Dq'^p^-'' to, respectively, i^)?^""^, {i)qp'~'^, and where 



1 



o + X 



-2[g,6]Re(^"') + {<?,o}Im(A'' 



-X^qdqB" 



(40) 



N2" = 1 + 2AgIm(A'") + X^q'^B" 



If qdq = 0, the terms {g2, o}Re(C""') and 6]Im(C"") 
should b e retai ned. If by chance one has also that 
{g2, 6} = i[q^,6\ — 0, one should in principle retain terms 
up to A'^, or higher, in the numerator. However, since 
by hypothesis A is small compared to the other relevant 
scales, one can neglect altogether these fine corrections, 
which lead to a small deviation of (o) from zero for ex- 
actly orthogonal pre- and post-selected states. 

Finally, for NOPPS, when AV^"" > AglA""], the 
leading term in the statistics is 



qe 



q- 



(41) 



This is a universal value, independent of the observable 
A. 



Zo{x\f)=N-'y2^-^ B. Large Ag* 



=0 



. / n\ . „ . When the term Xq* is large, with q* identifying a max- 

{-iy( . (r-^e^^°g^)o<„_„ (37) 

imum of the initial probe distribution \pQ{q,q')\ (typi- 
-'=0 cally, the initial average value q should be comparable to 



n 

E 
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- Intapolatit>n Simpler Inteipolationl 



(a) 



P. 

X 




Q 



• Exact Interpdati&n Simpler Inteipolationl 

(b) 



FIG. 2: Comparison between the interpolating formula 



Eq. (39 1, the simpler Eq. (401, and the exact result for the 



output p of a probe weakly measuring a spin component n er, 
as a function of the angle between the preselected state \ni^) 
and the postselected one | n/if) , n/ • rii = cos for a fixed pre- 
selection, with n ■ rii =0 and all vectors coplanar. The figures 
correspond to different interaction strengths: (a) \/5p — 0.5, 
(b) A/5p = 0.1. 



q* ) , the formulas given in the previous three sections the 
weak values must be substituted virith the expressions of 
Eqs. (|12b| and (flSbl), while q ^ q - q* . 



In particular, for an initial gaussian state of the detec- 
tor, q* = q, so that the even moments of p are 



(pi) f ^p3 + 



A2 (^2)p^~M5"' +Re(C"")] 



2 1 + A2(g - g)2 [B^ - Re(C»)] ' 



(42) 



and the odd ones 



{p>)f =Xj 



l-KA2(g-g)2[B«'-Re(C"' 



(43) 



VII. EXAMPLE: APPLICATION TO A SPIN 1/2 
OR Q-BIT 

We apply the results to the weak measurement of a 
spin component A = n- cr, for which there is an analytic 
solution [3U]. The exact conditional probability for the 
readout variable p is 

^(^1/) ^Z^ E K + <i + 2aRe«)] PoiP -\<y,p- \<J) 



[a^ - a^i + 2zcrlm(a5")] po(p + Xa,p- Act), 

(44) 



the normalization being 



N ■ 



1 + cos (2A(?) 



sin (2A(7)Im(a5" 



1 — cos (2Ag) 



The normal weak values are 

,„ 1 
2 
1 



a, 







(1 + n, • n/) 



Off =-n • (rij + n/ + x n/) , 
=^ (1 - n» • n/ + 2n • n; n • n/) , 



(45) 

(46a) 
(46b) 
(46c) 



where , are the polarizations of the pre- and post- 
selected density matrix, i.e., pij = (1/2)[1 -|- n^j • cr]. 
Equations (44) and (45) are new general results, since so 



far in the literature only an initial Gaussian state was 
considered for the probe [30l [S^ . 

For definiteness, we take the probe to be ini- 
tially described by the Wigner function Wo{p,q) oc 
exp{-p2/2Ap2 _ 2Jp2(^ _ q)2|^ where Sp < AP rep- 
resents its coherence scale. Then, the exact probability 
distribution for the rescaled readout A = p/\ is 



nA\f) 



1 



J2\ K+<i+2aRe«)] e" 
+ [cos (2Ag) {a^ - -f 2sin (2Ag)Im«)] 

(47) 



X e 



with A = AP/A, S = Sp/X, and 
1 + cos (2Ag)e-i/2'5' 



N ■ 



a^' + sin (2Aq)e-i/25 iin(a5") 



1 - cos(2Ag)e-i/2'5' 



-a 



(48) 



^ 1,1- 
The exact moments are, for even j, 

k=l ^ ^ 
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FIG. 3: Comparison between the interpolating formula (solid 
symbols) Eq. (22 1, and the exact result (unfilled symbols) 
Eq. (491 for the output (p")/ — p" in units of p" as a func- 
tion of the angle between the preselected state In^t) and the 
postselected one |n/:t), n/ • rii — cos 9 for a fixed preselec- 
tion, with n • Hi =0 and all vectors coplanar. Interaction 
strength: X/Sp = A/AP = 0.1. For the values of 6 not shown 
the discrepancy is undiscernible to the naked eye. 



and, for odd j, 



FIG. 4: Universal scaling for the moments. 



for even j, and, for odd j. 



{p')f -^Xjpr 



:j-Re(a5" 

N2 



(52) 



As Figs. [2] and [3] show, the interpolation formulas work 
satisfactorily over the whole range of pre- and postselec- 
tion. We notice that the validity of the interpolation is 
especially sensitive to the value of \q, since it contributes 
to first order. 

We notice that the even moments, to lowest order in 
A, present a universal scaling 



4A2Ar 



(53) 



Since for NOPPS af_i > and N ~ {l/25Ya^^^, all 
renormalized even moments take the value 6p^/AP'^ for 
orthogonal pre- and post-selection, as shown in Fig. |4j 
where we compare j = 2, 100, 400, showing that the scal- 
ing breaks down at n* ~ (AP/A)^. Odd moments, in- 
stead, vanish for orthogonal pre- and post-selection, since 
the second order in A in the numerator is zero. In prin- 
ciple, one should retain the third, or higher order terms, 
but this would give a value o{X/AP) for {p-') f, hence the 
interpolation is still satisfactory. 



TV 



0-l)/2 

E 

k=0 



(50) 



For even j, the interpolation formulas are given by 
Eqs. (22) and (23), with the a™ given above, or, equiv- 
alently, by the first non vanishing order expansion of 
Eqs. (|49| and (ISOl) 



{p')f~pi+X' 



nj-2. 



'1,1 



2No 



(51) 



VIII. CONCLUSIONS 

In conclusion, we have showed that the statistics of 
the weak measurement is interpolated satisfactorily by 
the rational function f /g, where both / and g are second 
order polynomials in the interaction strength A. The co- 
efficients are written in terms of averages over the initial 
state of the probe, and of a complex number, the canon- 
ical weak value, = Tr {pfApi}/ Ti {pfpi}, plus a real 
number = Tr {pfApi A}/ Ti {pfPi}. In some cases, 
the coefficients of Re{A^) and lm{A^) may vanish in the 
numerator /, so that it is necessary to introduce another 



9 



complex number, C" = Tv {pfA^pi}/ Tv {pf pi}, in order 
to obtain a satisfactory interpolation. 

The results presented here apply to an instantaneous 
interaction, so that the Hamiltonian evolution of the de- 
tector can be neglected. If this is not the case, the coef- 



ficients of the expansion have more involved expressions, 
and will be considered in a forthcoming work [36]. 

This work was supported by Fundagao de Amparo a 
Pesquisa do Estado de Minas Gerais through Process No. 
APQ-02804-10. 
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